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Abstract.  The  evolution  of  three-dimensional  disturbances  in  an  incompressible  three- 
dimensional  stagnation-point  flow  in  an  inviscid  fluid  is  investigated.  Since  it  is  not  possible  to 
apply  classical  normal  mode  analysis  to  the  disturbance  equations  for  the  fully  three-dimensional 
stagnation-point  flow  to  obtain  solutions,  an  initial-value  problem  is  solved  instead.  The  evolution 
of  the  disturbances  provide  the  necessary  information  to  determine  stability  and  indeed  the  com¬ 
plete  transient  as  well.  It  is  found  that  when  considering  the  disturbance  energy,  the  planar 
stagnation-point  flow,  which  is  independent  of  one  of  the  transverse  coordinates,  represents  a  neu¬ 
trally  stable  flow  whereas  the  fully  three-dimensional  flow  is  either  stable  or  unstable,  depending 
on  whether  the  flow  is  away  from  or  towards  the  stagnation  point  in  the  transverse  direction  that  is 
neglected  in  the  planar  stagnation  point. 
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1.  Introduction. 


Previous  analyticaJ  work  investigating  the  stability  of  planar  stagnation-point  flows  has  con¬ 
centrated  on  the  mathematical  simplification  provided  by  classical  mode  analysis  of  slreamwise  dis¬ 
turbances.  Wilson  and  Gladwell  (1978)  have  shown  that  incompressible  planar  stagnation  flow  is 
always  stable  to  three-dimensional  normal  mode  self-similar  disturbances  that  decay  exponentially 
outside  of  the  viscous  boundary  layer.  Lyell  and  Huerre  (1985)  re-examined  the  planar  stagnation 
flow  problem  by  using  the  same  class  of  disturbances  and  verified  the  results  of  Wilson  and 
Gladwell  as  well  as  characterized  the  other  stable  eigenvalue  branches  by  showing  that,  after  the 
initial  branch  found  by  Wilson  and  Gladwell.  the  other  branches  come  in  pairs.  In  addition,  by  a 
nonlinear  analysis  using  a  Galerkin  method,  these  authors  indicated  that  this  same  flow  is  unstable 
for  disturbances  of  sufficiently  high  amplitude.  On  the  other  hand,  a  numerical  study  by  Spalart 
(1989)  found  no  such  instability  suggesting  that  this  was  not  the  case.  Brattkus  and  Davis  (1991) 
showed  that  the  normal  mode  self-similar  disturbances  were  the  least  stable  of  the  class  of  distur¬ 
bances  that  have  a  power  like  behaviour  in  the  downstream  coordinate. 

Lasseigne  and  Jackson  (1992)  allowed  for  density  variations  induced  by  a  temperature 
difference  between  the  freestream  and  the  plate  and  determined  that  the  stagnation  flow  remained 
stable  to  small  streamwise  disturbances  regardless  of  the  plate  temperature.  The  effect  of  cooling 
the  plate  was  to  decrease  the  decay  rate  (less  stable)  of  the  small  wavelength  disturbances  while 
increasing  the  decay  rate  (more  stable)  of  the  moderate  wavelengths.  Again,  only  three-dimensional 
normal  mode  self-similar  disturbances  that  decay  exponentially  outside  of  the  viscous  boundary 
layer  were  considered.  Studies  dealing  with  the  swept  attachment  line  also  concentrated  on  stream- 
wise  normal  mode  linear  disturbances.  Hall,  Malik  and  Poll  (1984)  determined  that  a  region  of 
instability  (in  frequency- wave  number  space)  associated  with  increasing  crossflow  exists.  In 
independent  investigations  Kazakov  (1991)  and  Lasseigne,  Jackson  and  Hu  (1992)  determined  the 
effects  of  surface  temperature  variations  on  this  region  of  instability;  the  latter  investigation  also 
allowed  for  the  effects  of  suction  or  blowing  at  the  surface. 

It  is  not  possible  to  analyze  the  stability  of  the  fully  three-dimensional  stagnation-point  flow 
using  normal  mode  streamwise  disturbances  since  the  disturbance  equations  do  not  admit  this  class 
of  disturbances  as  an  eigenvalue  problem.  Therefore,  the  approach  taken  in  this  investigation  is 
different  than  that  of  previous  investigations.  A  more  general  initial  value  problem  is  solved  using 
the  methodology  developed  in  Criminale  and  Drazin  (1990)  and  has  its  origins  from  the  work  of 
Kelvin  (1887)  and  On  (1907a,b).  The  disturbances  are  taken  to  be  initially  bounded  in  all  direc¬ 
tions  and  the  evolution  of  initial  conditions  is  determined  analytically  and  in  closed  form.  By  con¬ 
centrating  on  the  mean  flow  subject  to  disturbances  in  an  inviscid  fluid,  the  fully  three-dimensional 
stagnation-point  flow  can  be  solved  with  the  planar  stagnation-point  flow  as  a  special  case. 

The  method  of  analysis  utilizes  a  moving  coordinate  transformation  that  allows  for  easy 
integration  (in  time)  of  the  individual  vorticity  components.  Then,  the  double  Fourier  transform  in 
the  new  transverse  coordinates  is  used  to  reduce  the  mathematical  problem  to  the  solution  of  ordi¬ 
nary  differential  equations  in  which  time  appears  strictly  as  a  parameter.  Thus,  a  completely  analyt¬ 
ical  solution  is  found  to  the  initial  value  problem  describing  the  evolution  of  threp  '^•mensional  dis¬ 
turbances  in  three-dimensional  stagnation-point  flow  in  an  inviscid  fluid.  The  time  evolution  of  a 
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single  Fnurier  mode  in  which  the  disturbance  is  periodic  in  the  trarcsverse  directions  is  investigated 
in  detail  as  well  as  the  evolution  of  an  initially  localized  disturbance.  In  both  cases,  the  evolution 
of  the  total  energy  of  the  disturbance  is  used  to  illustrate  the  importance  of  the  transient.  Kelvin 
(1887)  and  Orr  (19()7a,b)  have  shown  for  different  flow  configurations  that  there  can  be  growth  in 
the  perturbation  energy  even  when  the  classical  mode  analysis  shows  that  the  flow  is  stable  for 
long  times. 

Farrell  (1989)  applied  this  approach  to  two-dimensional  disturbances  within  a  two- 
dimensional  counterflow  as  a  good  approximation  to  the  local  flow  in  regions  of  confluence  and 
diffluence.  It  was  determined  that  plane  wave  disturbances  with  dependence  in  tlie  transverse  coor¬ 
dinate  and  independent  of  the  other  coordinate  have  energy  that  grows  exponentially  in  time.  Plane 
wave  disturbances  with  dependence  m  both  coordinates  experience  an  initial  exponential  growth, 
but  eventually  decay  as  time  progresses.  Disturbances  with  finite  wave  trains  that  were  not  spatially 
symmetric  were  shown  to  be  stable  but  these  disturbances  also  experienced  an  initial  transient 
growth  in  energy  with  the  energy  asymptotically  approaching  a  constant  amplitude.  Symmetric 
finite  wave  train  disturbances  were  shown  to  not  experience  the  initial  transient  growth  and  to  have 
energy  constant  in  time.  This  was  seen  to  occur  in  both  regions  of  confluence  and  diffluence. 

The  governing  equations  for  the  three-dimensional  stagnation-point  flow,  the  moving  coordi¬ 
nate  transformation  and  the  linear  disturbance  equations  are  presented  with  the  disturbance  equa¬ 
tions  solved  by  the  use  of  Fourier  transforms  in  the  transverse  coordinates.  Selected  results  for  the 
evolution  of  a  single  Fourier  mode  are  given  in  Section  3  and,  in  Section  4,  results  pertaining  to 
the  time  evolution  of  a  finite  wave  train  described  by  an  initially  Gaussian  profile  are  presented.  In 
Section  5,  a  constant  pressure  boundary  condition  is  considered  as  an  alternative  to  the  zero  normal 
velocity  condition.  Section  6  contains  a  discussion  on  the  effects  of  background  rotation  and  parti¬ 
cle  paths.  Conclusions  are  given  in  Section  7. 

2.  Problem  Statement  and  Basic  Equations. 

The  problem  under  investigation  is  that  of  linearized  disturbances  in  a  three-dimensional 
stagnation-point  flow.  The  basic  flow  is  given  nondimensionally  by 

U=x,  V  W='Kz,  (1) 

where  X  is  a  measure  of  three-dimensionality.  For  X  =  0  the  flow  is  a  two-dimensional  stagnation- 
point  flow;  for  X  =  1  the  flow  is  axisymmetric;  for  -0.4294  <  X  <  0  the  flow  corresponds  to  two 
symmetrically  displaced  protuberances  (Davey,  1961)  displayed  in  Figure  1.  For  X  <  0,  it  is 
important  to  observe  that  the  flow  is  toward  the  stagnation  point  in  the  z  -direction  and  away  from 
the  stagnation  point  in  the  -direction.  In  addition,  Davey  has  shown  that  separation  occurs  at 
X  =  -0.4294,  reversed  flow  exists  for  - 1  <  X  <  -0.4294,  and  no  solutions  are  possible  for 
X  <  - 1. 

The  nondimcib'onal  linearized  equations  for  small  disturbances  are  written  as 

Uj^  Vy  +  =  0 


(2a) 
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U,  +  XUjc  -i\  +X,)y  My+X2Uj+U  = 


(2b) 


V,  +  X  Vj,  -  (1+X,)>  Vy  +  Xz  Vj  -  (1+X)  V  =  -Py 


(2c) 


w,  +  X  Wj,  -  (1  +  X)y  Wy  +  Xz  Wj  +  Xw  =  -pj 


(2d) 


where  u,v,w  and  p  are  the  velocity  and  pressure  perturbations,  respectively.  The  appropriate 
boundary  conditions  require  that  all  disturbance  quantities  vanish  as  y  ->  <»  and  the  normal  velo¬ 
city  V  is  zero  on  the  wall.  In  addition,  initial  conditions  consistent  with  the  boundary  conditions 
must  be  supplied. 


The  above  equations  can  be  recast  in  terms  of  the  vorticity  components  in  the 

X ,  y ,  z  directions  and  are 


Dt 


=  ©r 


D©y 

“dT 


-(1+X)©y  , 


D©J 

~5r 


=  X©J . 


(3) 


where  D  IDt  is  the  material  derivative  defined  by 

=  ( )*  +  ( )jt  -  (1  +  X)y  ( )y  +  X z  ( )j .  (4) 


In  general,  only  two  vorticity  components  can  be  specified  at  time  r  =  0;  the  third  component  is 
found  by  appealing  to  the  equation  of  conservation  of  vorticity,  given  by 


d(Or  3©„ 


3x 


dy 


3©j 

3z 


=  0. 


(5) 


In  this  study,  the  initial  profiles  for  the  vorticity  components  ©^  and  ©^  are  specified,  and  the 
above  equation  is  used  to  determine  ©y.  Once  the  vorticity  has  been  determined  from  (3),  the  v 
velocity  component  is  found  from  the  following  relationship 


V^v 


d©; 

3x 


d(Ox 

dz 


(6) 


where  is  the  three-dimensional  Laplace  operator.  The  other  velocity  components  are  detennined 
by  af5)ealing  to  the  vorticity  relations  together  with  the  continuity  equation,  yielding: 


dxdy 


(7) 


d^v 

dydz 


(8) 


where  V2  is  the  two-dimensional  Laplace  operator  in  the  x  -  z  plane. 
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By  following  Criminale  and  Drazin  (1990).  an  analytical  solution  is  readily  obtained  if  a 
change  from  an  Eulerian  to  a  moving  coordinate  system  is  made.  The  moving  coordinate  transfor¬ 
mation 

^  =  Ti  =y  C  =  T  =  t  (9) 


is  chosen  so  that  the  system  of  equations  (3)  have  coefficients  that  are  functions  of  time  only  and 
the  material  derivative  (4)  becomes 


DO  ^  di) 
Dt  dT  ' 


(10) 


The  notable  feature  of  this  transformation  is  that  the  partial  differential  equations  for  the  individual 
vorticity  components  (3)  can  then  be  immediately  integrated  in  time  and  the  solutions  are 


0  T 

(0^ 


(0, 


0), 


(11) 


where  the  superscript  0  denotes  the  specified  values  at  T  =  0.  Once  the  initial  values  (0°  and  (0° 
have  been  made,  the  initial  value  of  0)°  is  determined  from  the  conservation  of  vorticity  (equation 
(5)).  The  equation  for  the  v  velocity  component  in  the  moving  coordinates  is 


a; 


(12) 


where  A  is  the  Ltqtlace  operator  in  the  moving  coordinates,  namely 


A  =  e 


-2T 


3^2 


+  e 


2(i+X)r 


an' 


+  e 


-2\T 


32 

ac' 


(13) 


It  is  noted  that  the  time  dependence  in  equation  (12)  appears  as  a  parameter  only  and  hence  finding 
a  solution  to  (12)  is  essentially  a  spatial  problem.  The  equations  for  the  other  two  velocity  com¬ 
ponents  in  the  moving  coordinate  system  are 


AtM  =  <»-0  +  2X)T 


3(0® 

"aT 


-  e 


XT 


3^v 

a^an 


(14) 


where 


3(0® 

'W 


+  e 


■2XT 


(15) 


(16) 


is  the  two-dimensional  Laplace  operator  in  the  moving  system. 

It  is  now  assumed  that  Fourier  transforms  may  be  taken  in  the  ^  and  ^  directions,  which 
implies  that  the  disturbances  are  bounded  in  these  directions.  The  double  Fourier  transform  is 
defined  by 
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j2(a.Ti.Y.T)  =  (J7) 

—  oo 

etc.  for  V,  w,  d>  with  the  inversion  given  by 

oo 

M(tn,CT)  =  — ^ffi2(a.Ti.Y.7')e-‘<“^^^^»rfadY.  (18) 

(Inydi 

The  transforms  of  the  relevant  equations  are 

^2(ux)r  _  5^2 .  ^  a<aO^-(i-^)r  +  ^  (19) 

3>}|0 

=  I  adj°  +  /  Y<J>°.  (20) 

dr\ 

ti=±  (yaO^-U-2X)T_,.„^^XT  ^  (21) 

a2  '  >  ail 


v0  =  -^  -iy~e^  .  (22) 

with 

G?  =  +  "f  .  (23) 

The  equation  for  v  is  a  second  order  ordinary  differential  equation  in  q  with  time  dependent 
coefficients,  and  since  Ot^  is  just  the  Fourier  transform  of  A2,  the  velocity  components  d  and  are 
determined  directly  from  the  algebraic  relations  provided  by  transforming  equations  (14)  and  (15). 

The  evolution  of  a  single  Fourier  mode  can  be  studied  by  choosing  the  initial  conditions  for 
the  vorticity  components  0)°  and  ® °  to  be 

(coi’.to")  =  (fl^  ,n,  )e"'^“”^^^“^^5(Ti-yo)  (24) 

since  (x ,  y ,  z )  =  (^,  q,  0  at  T  =  0.  In  Fourier  space, 

( <I)° , (Df)  =  ( , n, ) 5(a-ao) 5(y- Yo) 5(q -yo) •  (25) 

The  initial  value  for  the  last  vorticity  component  is  found  by  transfoiming  equation  (19)  and 
integrating  in  the  q-direction: 

(by  =  (i  aSl;,  +  i  Y^z  )5(a  -  ao)5(Y-  Yo)  [  fi(y\  -  yo)  -  1  ] 


(26) 
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where  H (r\)  is  the  step  function.  The  constant  of  integration  has  been  chosen  so  tJial  all  of  the 
vek)city  components  vatiish  as  ri  — >  c«. 


With  the  vorticity  components  known,  the  diftcrential  equation  for  v  is  integrated,  and 


V’  = 


+  X)V 
2C( 


/  ai2,,  -  /  Yi2,  t 

I  n-  viil  _ 


J 


5(a  -  ot(,)  g; ;  -  v„) 


-ft. 


tT1+>„) 


(27) 


where  &  is  as  defined  earlier.  The  solution  decays  as  r)  — ^  c«  and  is  zero  at  the  wall.  Once  the 
transform  is  inverted,  the  velocity  component  is  given  by 


2(2jt)‘Ct(,  I  J 


-ft„f 


(1  ♦Xifi 


n-.vol  _ 


(2X) 


where 


•  -IT  —'y'kT 

C(<,  =  UoC  +  Yo<’ 


(29) 


The  u  and  w  components  are  detennined  directly  from  the  earlier  equations  and,  upon  inverting 
the  transforms  these  components  are  found  to  be 


-Yo 


(27t)2&2  L 
«o 

2(271)^2 


00^2.*  +  Yo^ir 


^-(i+2X)r  ^-'(oo^+ro?) 


H(r]-yo)-l 


-yoQ^e 


-(I  +  2X)r  I  ^-'(Oo^  +  ToO 


(30) 


and 


w 


(^,T1,  T)  =  [  00^2^  +  Yo«r  ]  e  '  J  //(T|  -  y^,)  _  1  j 


Yo 


2(2jt)2cto^ 


oiijSl^e  -  Yo£2^  e 


[  ^-'(«o4+ro?) 


(31) 


3.  Single  Mode:  Infinite  Wave  Train. 

In  this  section  we  present  selected  results  for  the  single  mode  solutions  given  above  for  « ,  v , 
and  w .  We  first  note  that  the  solutions  are  a  linear  combination  of  £2^  and  £2^ .  Therefore,  there 
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exists  two  modes,  defined  as  follows: 

Mode  I  :  fl,  =  0.  ft,  =  1 

Mode  II  :  tJj  =  0 


The  results  for  Mode  I  are  shown  in  Figure  2  for  three  values  of  the  wavenumber  in  the  :  - 
direction  (Yq  =  0,  Yo  =  01-  Yo  =  10)  ^  single  value  of  the  wavenumber  in  the  jc -direction 

(Oo  =  10).  Plotted  in  this  figure  is  the  maximum  amplitude  of  the  velocity  components  I  u  I 
I  V  I  and  I  w  I  „„  as  a  function  of  time.  The  maximum  in  the  m  and  w  velocities  are  at  the 
wall,  while  the  maximum  in  the  v  velocity  is  in  the  interior  for  all  times  since  it  is  required  (hat  v 
vanish  at  the  wall.  In  each  sub-figure,  four  values  of  the  three-dimensionality  constant  were 
chosen:  a)  axisymmctric  flow,  X  =  1  (shown  in  the  figures  as  a  small-dashed  curve),  b)  X  =  0.5.  c) 
planar  stagnation-point  flow,  X  =  0  (shown  in  the  figures  as  a  large-dashed  curve),  and  d) 
X  =  -0.2. 

Included  in  these  figures  are  the  results  for  flie  special  cases  of  two-dimensional  disturbances 
in  a  two-dimensional  flow  (Yo  =  0,  X  =  0),  two-dimensional  disturbances  in  a  three-dimensional 
flow  (Yo  =  0,  X  ^  0),  and  three-dimensional  disturbances  in  a  two-dimensional  flow  (Yo  ^  0,  X  =  0). 

As  seen  from  the  figures,  all  three  components  of  the  flow  exhibit  exponential  decay  in  the 
long  time  solution  for  the  response  of  this  mode.  The  results  for  the  u  velocity  show  that  the  decay 
rate  is  independent  of  the  variable  X  and  of  the  wave  number  in  the  z  direction.  Since  the  analytic 
solution  for  the  velocities  is  known,  this  information  can  be  found  by  determining  the  behaviour  of 
equation  (30)  as  T  ->  «  and  is 

u  (^,  0,  C,  7-)  =  — ^  (32) 


where  the  first  term  of  (30)  is  dominant  when  X  <  1,  and  there  is  a  balance  between  the  first  and 
second  terms  for  the  axisymmetric  case  X  =  1  that  does  not,  however,  change  the  limit.  For  the 
long  time  behaviour  of  the  v  component  of  the  velocity  the  maximum  must  be  determined.  The 
maximum  is  found  to  always  occurs  at  q  =  yo  so  that 


C.7’) 


^  “oYo  -(3+X)r  («o^  +  7o?) 

(2tc)2 


(33) 


with  the  axisymmetric  case  again  being  a  special  limit  that  does  affect  the  amplitude  in  the  above 
equation  but  not  the  decay  rate.  The  asymptotic  behaviour  for  the  w  component  of  the  velocity  is 
also  determined  analytically.  From  equation  (31)  it  is  found  that 

H-(to.c.r)  =  -  -*^4--  - (34) 

(2Jt)^lYol 


The  dependence  of  the  decay  rate  of  the  v  and  w  components  on  the  three-dimensionality  parame¬ 
ter  X  is  clearly  seen  in  Figure  2. 
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Although  the  above  analytic  tbnns  appear  to  be  relatively  simple,  the  use  of  the  traiisfomicd 
variables  ^  and  p  tend  to  hide  some  rather  significant  changes  in  tlie  spatial  structure  of  the  dis¬ 
turbances.  For  X>  0  there  is  a  stretching  of  the  initial  Fourier  mode  fomi  in  both  the  v -direction 
and  the  r -direction.  For  the  axi.symmetric  ca.se,  ’his  .stretching  occurs  at  an  equal  rate  in  both  direc¬ 
tions.  For  the  planar  stagnation-point  flow  A.  =  0  there  is  lu)  stretching  of  the  disturbance  in  the  r 
direction,  and  for  A.  <  0  where  the  mean  flow  is  towards  the  stagnation  point  in  the  z -direction, 
there  is  an  actual  contraction  of  the  irtitial  Fourier  mode  form  in  the  r  -direction.  In  the  y  direction, 
it  is  noted  that  although  the  maximum  in  v  velocity  is  at  a  fixed  value  of  r).  this  implies  mat  this 
maximum  approaches  the  wall  expoirentially  fast  when  the  problcnr  is  converted  to  the  physical 
spatial  variables,  that  is  the  inverse  of  (9). 

Although  the  long  time  behaviour  of  the  disturbances  is  id  course  very  important,  equally 
important  is  the  transient  evolution  since  it  is  possible  for  there  to  be  signilicant  initial  growili  (or 
perhaps  a  lotig  time  persistence)  of  the  disturbance  quantities  before  the  inevitable  exponential 
decay.  By  examining  the  results  of  Figure  2,  it  is  seen  that  the  initial  response  to  a  Mode  I  di.stur- 
bance  is  typically  a  linear  decay  in  time.  The  initial  linear  decay  rate  depends  on  the  value  of  the 
three-dimensionality  parameter  X  with  X  <  0  decaying  at  the  slowest  rate.  The  X  dependence  of  this 
initial  decay  rate  is  strongest  for  smaller  wave  numbers  ir  the  z -direction. 

The  results  for  the  single  mode  response  to  a  Mode  II  disturbance  are  shown  in  Figure  A 
significant  difference  from  the  Mode  I  results  is  immediately  noticeable.  Analytically  the  long-time 
limits  are  given  by 

„(5.0,c,r)  ,-2(i.X)r  ,-.  (.^1.7.?) 

(2  n)'  I  Yo  I 

and 

(2  It)' 

as  7  — >  oo.  It  is  .seen  from  Figure  3  that  the  u  and  v  components  of  velocity  decay  exponentially 
for  all  X  while  the  w  component  of  velocity  only  decays  when  the  three-dimensionality  parameter 
X  is  greater  than  zero.  This  value  of  X  con  -londs  to  a  stretching  of  the  original  disturbance  in  the 
z  direction.  However,  for  the  planar  stagnation-point  flow  where  X  =  0  the  w  component 
approaches  a  constant  independent  of  the  initial  wave  number  in  the  x  and  z  directions.  For  the 
case  where  there  is  mean  flow  towards  the  stagnation  point  along  the  z -direction  (X  <  0),  the  w 
component  grows  exponentially.  It  is  seen  that  the  initial  behaviour  of  the  u  and  v  components 
also  reflect  this  differing  behaviour  for  the  planar  stagnation-point  flow  and  the  X  <  0  case.  It  is 
also  true  that  there  is  initially  linear  growth  of  these  components.  As  in  the  Mode  I  case,  the 
dependence  of  the  behaviour  for  ear.y  time  on  the  three-dimensionality  parameter  is  strongest  at 
smaller  wave  numbers  in  the  z -direction 


(35) 

(36) 

(37) 
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The  results  presented  in  Figures  2  and  3  give  the  maximum  values  of  the  velocities,  but  con¬ 
sidering  the  severe  contraction  of  the  disturbances  in  the  y  -direction  as  evidenced  by  the  movement 
of  the  maximum  in  the  v  velocity  toward  the  wall  at  an  exponentially  fxst  rate,  it  is  perhaps  nu)rc 
useful  to  kH)k  at  the  growth  or  decay  of  the  energy  of  the  disturbances.  For  the  single  Foimer 
nuKle  results  presented  in  this  section,  it  is  necessary  to  consider  the  energy  per  period  in  the  c,  and 
s  directions,  or  equivalently  the  energy  per  period  in  the  v  and  r  directions.  Because  of  the 
stretching  in  the  .v -direction  and  the  stretching  or  contraction  (depending  on  the  sign  of  X)  iti  the 
: -direction,  these  results  do  not  readily  lend  themselves  to  direct  physical  interpretation;  however, 
they  do  help  point  out  the  effect  that  the  contraction  has  in  the  vertical  direction  v .  Shown  in  Fig¬ 
ure  4  is  the  quantity  £  as  a  function  of  time  defined  by 


d\ 


(38) 


where  u.  v.  and  w  are  given  by  (28),  (30),  and  (31).  Each  graph  was  normalized  by  the  value  at 
time  t  =  0.  The  amplitudes  and  i2,  are  equal  in  the  mixed  mode  results.  In  each  sub-figure, 
four  values  of  the  three-dimensionality  constant  were  chosen:  a)  axisymmetric  flow.  X  =  1  (shown 
in  the  figures  as  a  small-dashed  curve),  b)  planar  stagnation-point  flow,  X  =  0,  c)  X  =  -  1  /  3 
(shown  in  the  figures  as  a  large-dashed  curve),  and  d)  X  =  -0.4.  From  these  figures.  Mode  11  is 
seen  to  represent  disturbances  that  may  grow  in  energy  as  it  also  represented  disturbances  with 
growth  in  the  maximum  values  of  the  velocities.  However,  it  is  noted  that  although  the  maximum 
in  the  velocities  can  grow  for  any  X  <  0,  the  energy  per  period  grows  only  for  X  <  -  1  /  3.  This 
phenomenon  is  directly  related  to  the  contraction  in  the  y  direction.  Very  little  dependence  on  die 
wave  number  in  the  z  direction  can  be  detected  which  is  consistent  with  the  results  of  Figures  2 
and  3.  The  energy  of  Mode  1  is  dominated  by  the  u  and  v  velocities  which  show  little  depen¬ 
dence  on  Yo-  ^d  the  Mode  II  energy  is  dominated  by  the  w  velocity  which  also  shows  little  depen¬ 
dence  on  Yo- 

4.  Finite  Wave  Packet. 


In  the  previous  section,  the  results  for  an  initial  disturbance  that  consisted  of  a  single  Fourier 
mode  were  presented.  This  disturbance  has  an  infinite  spatial  profile  in  the  x  -  z  plane.  Considering 
the  distortions  in  the  x  and  z  directions  introduced  by  the  transformed  variables,  there  is  consider¬ 
able  doubt  as  to  the  proper  interpretation  of  the  single  mode  results.  To  resolve  this  difficulty  in 
the  analysis,  an  initial  disturbance  that  is  localized  in  space  is  chosen.  The  initial  vorticily  profile  is 
taken  as  a  Gaussian.  This  disturbance  is  initially  symmetric  in  the  x  -  z  plane  and,  since  the  invis- 
cid  problem  is  being  investigated,  a  Gaussian  profile  in  the  y  direction  is  also  chosen  and  all  boun¬ 
dary  conditions  can  still  be  satisfied. 


The  analysis  proceeds  by  replacing  (24)  with 


(0)“,(0°)  = 


_2li 

4P 


e 


(39) 


-  lo¬ 


in  Fourier  space. 


_Jii 


(<1,0  ^0^^  a(a^+Tf=)  4P 

(27rVjir 


(40) 


where  the  constant  is  chosen  such  that  the  integral  of  the  square  of  the  vorticity  over  the  domain  is 
a  constant  for  each  choice  of  o  and  p. 

The  above  initial  conditioas  can  be  substituted  directly  into  equations  (19)  and  (20)  which 
can  then  be  solved  for  v  and  or  the  previous  results  for  a  delta  function  initial  condition  in  the 
y  direction  can  be  integrated  over  the  dummy  variable  Vo-  Either  way,  the  transform  of  the  vorti¬ 
city  0)^’  and  the  v  component  of  velocity  are  given  by 


<  =  -(/■  a  LI,  +  i  yil- ) 

■  (27tVp)''^ 


2Vp 


(41) 


and 
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_  -  (1  ♦  Mr 


n_  +  a^-(ux.7-vp 
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Cte 


-o+mt 
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.  (42) 


The  transforms  of  the  u  and  w  components  are  found  directly  from  equations  (21)  and  (22)  and  all 
transforms  can  be  inverted  using  the  formula  (18).  The  long  time  behaviour  of  w  at  the  wall  for 
Mode  II  is  given  by 


(27cVp)‘"‘ 


(43) 


as  T  which  shows  that  the  transform  grows  exponentially  for  X  negative  consistent  with  the 
results  of  the  single  mode. 

In  order  to  exantine  the  true  behaviour  of  the  evolution  to  this  initially  localized  disturbance, 
the  energy  as  a  function  of  time  is  computed.  The  energy  in  terms  of  the  unstretched  variables  is 


1 


£■(7)  =  -j  I  JJ  [u^  +  v^  +  yv^]dx  dz  dy. 


(44) 


which  can  also  be  given  in  terms  of  integrals  over  the  Fourier  transform  quantities  in  the  stretched 
variables  as 


EiT)  = 


,(i+X)r 
8  71^ 


,ye^^ ,y  ,T)dU.dydy , 


(45) 
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where 

Q\a,y,r\,T)  =  u^+v~  +  w".  (46) 

The  normalized  energy  of  the  response  to  a  Mode  II  gaussian  profile  with  P  =  1.0  and  a  =  0.5  is 
shown  in  Figure  5  for  eight  values  of  the  three-dimensionality  parameter  X  ranging  from  X  =  1  (the 
axisymmetric  case)  to  X  =  -  .4  (near  separation).  A  number  of  interesting  features  are  noticeable. 
First,  for  the  planar  stagnation- point  flow  (X  =  0).  the  energy  approaches  a  constant  value  approxi¬ 
mately  three  times  the  initial  value.  This  is  in  agreement  with  the  results  tracking  the  maximum  in 
the  velocities.  The  mean  flow  produces  a  contraction  of  the  disturbance  in  the  y  -direction,  initially 
increasing  the  magnitude  of  the  vv  velocity  at  the  wall  and  initially  increasing  the  energy  of  the 
disturbance.  However,  a  balance  develops  between  the  continued  contraction  in  the  y -direction  and 
the  relief  provided  to  the  flow  in  the  .v -direction  thereby  leading  to  the  situation  where  the  energy 
approaches  a  constant.  This  behaviour  of  the  disturbance  ettergy  for  the  two-dimensional  planar 
counter-flow  was  also  seen  by  Farrell  (1989).  The  purpose  of  this  study  is  to  determine  the  effects 
of  the  three-dimensionality  of  the  mean  flow.  When  the.se  three-dimensional  effects  are  considered 
(X  *  0),  it  is  seen  that  the  planar  stagnation  point  flow  represents  a  special  case.  For  X  slightly 
positive,  there  is  again  an  initial  transient  growth  in  energy.  The  extra  relief  provided  to  the  flow 
in  the  r -direction  means  that  the  balance  that  developed  for  the  planar  stagnation-point  flow  does 
not  develop  and  the  energy  of  the  disturbance  eventually  decays  to  zero.  For  X  sufficiently  large 
and  positive  there  is  no  initial  transient  growth  implying  that  the  relief  provided  to  the  flow  in  the 
r -direction  is  sufficient  to  prevent  the  the  initial  increase  in  disturbance  energy.  From  Figure  5,  it 
is  .seen  that  for  the  .special  case  of  axisymmetric  flow  (X  =  1),  the  energy  undergoes  an  immediate 
exponential  decay.  For  X  less  than  zero,  an  unstable  situation  develops,  and  the  energy  of  the  dis¬ 
turbance  continually  increases  at  an  exponential  rate  owing  to  the  contraction  in  both  the  y  and  : 
directions  which  is  not  balanced  by  the  relief  in  the  x  direction. 

It  is  also  interesting  to  examine  the  behaviour  of  the  vorticity  components  which  are  given  by 
equation  (11)  in  the  transformed  coordinates.  In  the  physical  coordinates,  the  vorticity  field  is 
undergoing  a  contraction  in  the  y -direction,  an  expansion  in  the  jr -direction  and  an  expansion  (for 
X  >  0)  or  a  contraction  (for  X  <  0)  in  the  2 -direction.  The  exponential  time  factors  in  equation  (11) 
indicate  that  for  a  Mode  II  disturbance  vorticity  is  transferred  from  the  y -component  to  the  .r- 
component  while  the  transfer  of  vorticity  between  the  y  and  2  components  for  a  Mode  1  distur¬ 
bance  depends  on  the  sign  of  X.  By  considering  the  response  of  a  Mode  II  initial  disturbance,  it  is 
seen  that  in  the  neutrally  stable  case  (X  =  0)  the  decay  rate  of  the  y -vorticity  component  is  equal  to 
the  growth  rate  of  the  .v  -vorticity  component.  For  the  stable  ca.se  (X  >  0),  the  decay  rate  of  the  y  - 
vorticity  component  is  greater  titan  the  growth  rate  of  the  ,v -vorticity  component;  whereas,  for  the 
unstable  case  (X  <  0),  the  decay  rate  of  the  y -vorticity  component  is  less  than  the  growth  rate  of 
the  X -vorticity  component.  For  the  Mode  I  disturbance,  either  both  components  are  decaying 
(X  <  0),  or  the  decay  rate  of  the  y  -vorticity  component  is  always  greater  than  the  growth  rate  of 
the  2 -vorticity  compt)ncnt.  Perhaps  this  is  an  indication  of  why  it  is  the  Mode  II  disturbances  that 
experience  uastable  energy  growth  rather  than  those  of  Mode  1. 
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5.  Constant  Pressure  Boundary  Condition. 

One  problem  of  interest  in  geophysical  fluid  dynamics  is  the  case  where  there  is  a  constant 
pressure  surface  in  lieu  of  a  condition  on  the  velocity.  Such  a  problem  has  been  examined  by  Eady 
(1949),  for  example,  where  plane  Couette  flow  is  mathematically  equivalent  when  the  constant 
pressure  surface  is  used.  Consider  p  =  0  at  y  =  0  rather  than  the  v  velocity.  Of  course  the  surface 
V  =  0  can  no  longer  be  thought  of  as  a  solid  wall.  To  derive  an  equation  which  governs  the  time 
evolution  of  the  pressure,  the  momentum  equations  (2b)  and  (2c)  can  be  combined  with  the  con¬ 
tinuity  equation  (2a)  in  the  (i^,r|,^,  T)  coordinate  system  to  get 

■■  oT  dT| 


where  it  is  immediately  seen  that  the  axisymmetric  case  A.  =  1  is  a  special  case.  In  keeping  with 
the  previous  analysis,  the  pressure  equation  is  first  transformed  into  Fourier  space  and  then  u  is 
eliminated  by  using  (21).  yielding 


=_  lia  (X_t)(i)%-2"+^'^ 


+ 


a* 


dv 

dq 


dr  ^  dn 


(48) 


If  the  initial  conditions  are  considered  to  be  given  by  equation  (24)  for  this  problem  so  that  only  a 
single  Fourier  mode  is  investigated,  then  fit®  is  given  by  (26),  while  (19)  can  be  integrated  to  find 
V,  yielding 


V  =  A  i(T)e 


-(i-X)ri 


Tl-yol 


+  A2{T)e 


-ttf 


(n->'o) 


(49) 


where 


Ai{T)  = 


.ll+X)?" 

2a 


/■  aiA,  \  5(a  -  cto)  5(y  -  Yq)  . 


(50) 


The  determination  of  A  i(7)  satisfies  all  of  the  conditions  at  q  =  yo  and  it  is  left  to  impose  the  con¬ 
dition  on  the  pressure  at  q  =  0  in  order  to  determine  the  remaining  unknown  function  Ai{T). 
Equation  (48)  provides  the  necessary  information  in  determining  A  2,  but  it  is  convenient  to  work 
with  an  alternate  dependent  variable  defined  by 


B{T)  s  -^(O.T)  = 

dq 


(i-Mr 


-  4-)t? 


ftvor 


(1  ♦A.ir 


(51) 


The  equation  for  B(T)  is  then 


ciT  or 


(52) 


where 
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ii  =  (/  aii,  +  /  )5(a  -  a,i)5( y-  Yo). 


(53) 


It  is  helpful  find  the  long  time  solution  in  order  to  detennine  the  effects  of  the  three- 
dime’'  »)nulity  parameter  X.  Since  Mode  II  (ii,  =  l.ii-  =  0)  is  the  more  interesting  initial  condi- 
ti<'  the  rest  of  the  analysis  is  restricted  to  this  case  with  X<  1.  The  long  time  solution  for  B  is 
determined  directly  fn)m  (52)  or 

B  =B„U-  (54) 

r 

as  T  oo  where  the  the  proportionality  to  i2  is  explicitly  given  and  the  coastant  B^  can  be  found 
numerically.  The  behaviour  of  A  ^  is  determined  directly  from  equation  (50),  and  A  2  from  equation 
(5 1 ).  These  limits  are 

Ai(T)  =  -^  5(a  -  a<08(y -  y„) .  (55) 

A2iT)=^--^i^ -aBJe'^'^  6ia-a^)b(y-y^,).  (56) 

I  y  I  2 

i\s  T  00.  The  long  time  behaviour  of  the  velocities  u,v,w  is  found  through  equations  (24). 
(49)  and  (22)  indicating 

M(0.r)=  4(Y-afi.)c-^6{a-ao)6(y-Yo).  (57) 

r 

V  (0. T)  =  -  -p^  (y  -  5(a  -  Oo) 6(y  -  yo) ,  (58) 

vv(0,r)  =  --^B„e-^^6{a-(Xo)6(y-yo).  (59) 

y 

as  7  00.  It  is  seen  that,  for  the  boundary  condition  p  =  0  on  the  surface  y  =  0,  the  velocity 

components  v  and  vv  of  Mode  II  are  exponentially  growing  when  the  three-dimensionality  parame¬ 
ter  X  is  negative. 

6.  Effects  of  Background  Rotation;  Particle  Paths. 

The  expression  used  for  the  basic  velocity  and  given  by  (1)  can  be  modified  to  include  effects 
due  to  backgrour.J  rotation  (strain).  Specifically,  if  is  the  constant  dimensional  rate  of  rotation, 
then 

U=x+Q„z,  V  =-(l-l-X,)>-,  W=A.z-n^.r.  (60) 

becomes  the  new  non-dimensional  representation.  The  process  of  shifting  to  a  moving  coordinate 
system  and  solving  for  the  perturbations  as  previously  done  can  likewise  be  accomplished  here  but 
with  a  noticeable  increase  in  complexity.  For  example,  the  fundamental  vorticity  equations  (3)  used 
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for  solving  the  initial-value  problem  are  no  longer  uncoupled  and  the  net  effect  is  a  higher  order 
differential  system.  Thus,  the  dynamics  is  altered  to  some  degree.  In  particular,  temporal  oscilla¬ 
tions  become  possible  but  the  overall  stability  conclusions  remain.  These  points  can  be  illustrated 
be  examining  the  paths  of  the  material  particles  in  the  basic  flow. 

The  Cartesian  components  for  any  particle  in  a  Lagrangian  frame  can  be  written  using  (60)  as 

—  =  U=x+il.z.  (61a) 

dT 

^  =  V=-(l+X)y,  (61b) 

dT 

—  =  w  =Xz  —  X  .  (61c) 

dT 


By  assuming  x  =  ,  y  =  and  z  =  at  time  T  =  0,  the  solutions  for  (61a,b,c)  are 

(62a) 


-(I+X)T 


(62b) 
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.XT 


oX 


(62c) 


with 


a  = 


1  +  X 
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X-t-Q^  ' 

1/2' 

1  -1- 

1-4 - ^ 

(1  +  > 

(63) 


Oscillatory  solutions  are  possible  if  Qo  >  y  (I  -  ^)- 
from  any  initial  position. 


In  general,  particles  will  move  arbitrarily  far 


Two  interesting  limits  of  (62a,b,c)  are  when  =  0  or  X  =  1;  X,  can  never  be  very  negative 
and  therefore  a  change  in  the  system  behaviour  using  this  parameter  is  not  possible.  In  the  first 
instance,  =  0,  then 

x=x„e^,  y  ,  z  =  z„  .  (65) 


The  novel  feature  of  (65)  is  that  these  results  are  identical  to  those  of  (II)  for  the  vorticity  where 
jc,y,z  are  replaced  by  (Oj^,  to^,,  (o^.  In  short,  this  is  a  unique  situation  where,  for  a  three- 
dimensional  flow,  the  particle  paths  are  synonymous  with  vortex  lines.  Extending  this  argument  to 
the  ^  0  case  can  be  made  by  conjecture  since  the  vorticity  has  not  bf^en  determined  under  these 
circumstances.  There  is  little  likelihood  that  the  unstable  aspects  of  the  flow,  however,  can  be 
changed  by  a  finite  When  X  =  1,  the  basic  flow  is  axisymmetric  and  a  =  1  -i-  /  indicating 
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that  background  rotation  has  caused  the  system  to  have  inertia  and  therefore  oscillations  in  time. 

7.  Conclusions. 

We  have  investigated  the  evolution  of  three-dimensional  disturbances  in  a  fully  three- 
dimensional  stagnation-point  flow  in  an  inviscid  fluid.  It  has  been  shown  that  the  planar 
stagnation-point  flow  is  a  special  case  in  which  the  disturbance  energy  appn)aches  a  constant  for 
long  time.  If  the  flow  in  the  second  transverse  coordinate  is  away  fmm  the  stagnation  point  then 
the  flow  is  provided  enough  relief  such  that  the  disturbance  energy  decays  after  an  initial  transient 
growth.  As  the  limiting  case  of  axisymmetric  flow  is  approached,  the  disturbance  energy  is  found 
to  decay  without  an  initial  transient  growth.  For  flow  towards  the  stagnation  point  in  the  second 
trairsverse  direction,  it  is  found  that  the  disturbance  energy  may  grow  exponentially  thus  indicating 
an  unstable  flow  configuration.  The.se  results  were  found  by  determining  a  closed  form  solution  to 
the  initial  value  prtrblem  even  though  a  cla.ssical  mode  analysis  was  not  possible  for  the  fully 
three-dimensional  flow. 

Because  of  the  inviscid  assumption,  the  results  for  the  planar  stagnation-point  flow  cannot  be 
compared  directly  with  previous  work  but  the  method  used  here  can  be  extended  (with  consider¬ 
ably  more  mathematical  complexity)  to  the  study  of  the  inviscid  mean  flow  subject  to  viscous 
linear  disturbances.  These  results  can  be  compared  with  .some  of  the  previous  work.  However,  any 
results  for  this  ptxrblem  can  only  be  suggestive  in  view  of  the  fact  that  the  basic  flow  is  derived 
from  an  inviscid  basis. 

Background  rotation  or  strain  of  the  field  can  only  alter  the  dynamics  by  forcing  the  system 
to  have  temporal  oscillations  rather  than  prevent  instabilities.  Finally,  more  comprehensive  material 
particle  path  information  could  be  obtained  if  one  allows  for  the  velocity  field  to  include  the  pertur¬ 
bations  as  well  as  the  basic  flow.  Pursuit  of  this  information  will  involve  three-dimensional  coupled 
nonlinear,  non-autonomous  differentia'  eauations. 
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Figure  1.  Schematic  three-dimensional  stagnation  point  flow. 
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Figure  2.  Maximum  amplitude  of  the  velocity  components  as  a  function  of  time  for  Mode  1 
disturbances  for  three  values  of  the  z -wavenumber.  In  each  sub-figure,  four  values  of  the 
three-dimensionality  constant  chosen  are:  a)  axisymmetric  flow,  k  =  1  (shown  in  the  figures 
as  a  small-dashed  curve),  b)  A,  =  0.5,  c)  planar  stagnation-point  flow,  X  =  0  (shown  in  the 
figures  as  a  large-dashed  curve),  and  d)  X  =  -0.2. 
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Figure  3.  Maximum  amplitude  of  liie  velocity  components  as  a  function  of  lime  for  Mode  11 
disturbances  for  three  values  of  the  z -wavenumber.  In  each  sub-figure,  four  values  of  the 
three-dimensionality  constant  chosen  are:  a)  axisymmelric  flow,  X  =  1  (shown  in  llte  figures 
as  a  small-dashed  curve),  b)  X.  =  0.5,  c)  planar  stagnation-point  flow,  X.  =  0  (shown  in  the 
figures  as  a  large-dashed  curve),  and  d)  X  =  -0.2. 
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Figure  4.  Nonnalized  energy  £  as  a  function  of  time.  In  each  sub-figure,  four  values  of  the 
Iliree-diinensionality  constant  cliosen  arc:  a)  axisymmelric  How,  X  =  I  (shown  in  the  figures 
as  a  smail-dashed  curve),  b)  planar  stagnation-point  (low,  >.  =  0,  c)  X  =  -  I  /  3  (shown  in  the 
figures  as  a  large-dashed  curve),  and  d)  X  =  -0.4. 
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